The data for the models was derived for the initial dataset dsL:

numID<- 9022 # highest id value (max = 9022)
### Define the data that will populate the model
ds<- dsL %>%  # chose conditions to apply in creating dataset for modeling
  dplyr::filter(id < numID) %.% # 1:9022
  dplyr::filter(year %in% c(2000:2011)) %.% # 1997:2011
  dplyr::filter(sample %in% c(1)) %.% # 0-Oversample; 1-Cross-Sectional
  dplyr::filter(race %in% c(4)) %.% # 1-Black; 2-Hispanis; 3-Mixed; 4-White
  dplyr::filter(byear %in% c(1980:1984)) %.% # birth year 1980:1984
  dplyr::filter(ave(!is.na(attend), id, FUN = all)) %.% # only complete trajectories
  dplyr::mutate( # compute new variables
    age= year-byear, # definition of age to be used in the model    
    timec=year-2000, # metric of time is rounds of NSLY97 in years, centered at 2000
    timec2= timec^2, 
    timec3= timec^3,
    #         timec= age-16, # metric of time is bilogical age in years, centered at 16
    #         timec2= timec^2,
    #         timec3= timec^3,# 
    cohort=byear-1980) %.% # age difference, years younger (unit - 1 cohort away)
  dplyr::select( # assemble the dataset for modeling
    id, sample, race, byear,cohort, # Time Invariant variables
    year,
    age, timec,timec2,timec3, attend)  # Time Variant variables

m0F – Fixed only

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m0F

m1F

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m1F

m2F

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m2F

m3F

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m3F

m4F

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m4F

m5F

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m5F

m6F

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m6F

m7F

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m7F

mFa –

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mFa

mFb

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mFb

mFc

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mFc

mFf

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mFf

mFd

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mFd

mFe

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mFe

m1R1 – 1 Random

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m1R1

m2R1

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m2R1

m3R1

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m3R1

m4R1

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m4R1

m5R1

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m5R1

m6R1

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m6R1

m7R1

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m7R1

mR1a –

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR1a

mR1b

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR1b

mR1c

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR1c

mR1f

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR1f

mR1d

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR1d

mR1e

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR1e

m1R2 – 2 Random

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m1R2

m2R2

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m2R2

m3R2

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m3R2

m4R2

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m4R2

m5R2

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m5R2

m6R2

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m6R2

m7R2

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m7R2

mR2b –

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR2b

mR2c

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR2c

mR2f

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR2f

mR2d

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR2d

mR2e

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR2e

m2R3 – 3 Random

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m2R3

m3R3

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m3R3

m4R3

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m4R3

m5R3

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m5R3

m6R3

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m6R3

m7R3

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m7R3

mR3f –

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR3f

mR3d

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR3d

mR3e

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk mR3e

m3R4 – 4 Random

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m3R4

m4R4

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m4R4

m5R4

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m5R4

m6R4

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m6R4

m7R4

Called from: BuildLine(modelName = modelName)
debug: p <- ggplot2::ggplot(ds, aes(x = timec, y = attend, group = id))
debug: p <- p + geom_line(aes(x = timec, y = yHat), colour = "red", 
    alpha = indLineAl, size = indLineSz, na.rm = T)
debug: p <- p + geom_line(aes(y = yFE), fill = NA, na.rm = T)
debug: p <- p + scale_x_continuous(breaks = seq(1, 11, 2))
debug: p <- p + scale_y_continuous(breaks = seq(0, 8, 1))
debug: p <- p + coord_cartesian(xlim = c(-0.5, 11.5), ylim = c(0.5, 
    8.5))
debug: p <- p + themeLine
debug: p <- p + labs(title = "How often have you attended a worship service?", 
    x = "Years Since 2000", y = "Church attendance")
debug: return(p)

plot of chunk m7R4

Model specification

General case of LCM and model sequence are defined as:
general form